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Abstract 

A new mechanism for the evolution of primordial magnetic fields is described and 
analysed. The field evolution is followed from the time of its creation until the epoch 
of structure and galaxy formation. The mechanism takes into account the turbulent 
behaviour of the early universe plasma, whose properties determine strongly the 
evolution of the field configuration. A number of other related issues such as the case 
of an electroweak plasma are also considered. Finally, as an example, the mechanism 
is applied to specific models. 
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1. One of the most exciting astrophysical consequences of phase transitions in the 
early universe is the possible creation of primordial magnetic fields. The existence of a 
primordial magnetic field could have significant effect on various astrophysical processes. 
In particular it may be involved in the galaxy formation process fl]], or in the generation 
of the observed galactic magnetic fields. 

It is widely accepted that the galactic magnetic fields are generated through a galactic 
dynamo mechanism. Here, a weak seed field is exponentially amplified by the turbulent 
motion of ionized gas, which follows the differential rotation of the galaxy [Q], 0]. The 
currently observed magnetic field of the Milky Way and of nearby galaxies is of the order 
of a ixGauss. If the e-folding time is no more than the galactic rotation period, ~ lO'^yrs, 
then, considering the galactic age, ~ lO^^yrs, the seed field needed to produce a field 
of the observed value is about ~ lQ~^'^Gauss on a comoving scale of a protogalaxy 
(~ 100 kpc). Since the gravitational collapse of the protogalaxies enchances their frozen-in 
magnetic field by a factor of (pc/Po)^^^ ~ 10'^, where pc ~ lO"^'^^' cm~^ is the typical mass 
density of a galaxy and po — 2 x g cm~^ is the current cosmic mass density, this 

seed field corresponds to an rms field of the order of ~ over the comoving 

scale of ~ 1 Mpc. With the assumption that the rms field scales approximately as 
with the expansion of the universe, where a oc t^/'^ is the scale factor in the matter era, we 
find that the required rms value of the seed field at the time teq ~ lO^^sec of equal matter 
and radiation densities is ~ 10~^^GaMss x {tgc/teqY/'^ ~ 10~^°GaMss, where tgc ~ lO^^sec 
is the time of the gravitational collapse of the galaxies [Q]. Thus, an rms field of magnitude, 
B^l^g ~ IQ^'^^Gauss at t^q would be sufficient to seed the galactic dynamo and generate 
the observed galactic magnetic fields. 

Various attempts have been made to produce a primordial field in the early universe 
0. In most of the cases, though, the achieved field appeared to be too weak or incoherent 
to seed the galactic dynamo. In this letter, it is shown that the key issue, determining the 
rms value of a primordial magnetic field at tgg, is the evolution of the correlated domains 
of the field, i.e. the growth of the lengthscale over which the magnetic field is coherent. 
We develop a detailed mechanism for the evolution of the magnetic field configuration and 
we show that, in contrast to what is usually assumed, the correlation length, in general, 
grows faster than the scale factor a. This results in more coherent rms fields at the epoch 
of galaxy formation. 

2. In the existing literature the entire magnetic field configuration is taken to be 
comovingly frozen. As a result the rms magnetic field is thought to evolve as a~^, due to 
flux conservation. 

The fact that the early universe is, with great accuracy, a perfect conductor, ensures 
that magnetic flux is, indeed, conserved (in absense of dissipation) and, therefore, the 
magnetic field can be thought to be "frozen" into the plasma §], 0; over a certain scale. 
However, the simplistic assumption that the correlated domains of the field expand only 
due to the Hubble expansion does not take into account that the faster expanding causal 
correlations, through electromagnetic turbulence, could rearrange the field and correlate it 
on comoving scales larger than its initial correlated domains. This is because, when two 
initially uncorrelated neighbouring domains come into causal contact, the magnetic field 
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around the interface is expected to untangle and smooth, in order to avoid the creation 
of energetically unfavoured magnetic domain walls. In time the field inside both domains 
"aligns" itself and becomes coherent over the total volume. The velocity f , with which 
such a reorientation occurs, is determined by the plasma, which carries the field and has 
to reorientate its motion for that purpose.^ 

Thus, the evolution of the correlation length is given by, 

(1) 

where ^ is the correlation length of the magnetic field configuration, H is the Hubble 
parameter and v is the peculiar, bulk velocity, determined, in principle, by the state of the 
plasma. 

From (|I|) it is apparent that the correlated domains could grow faster than Hubble 
expansion. Therefore, the magnetic field configuration is not necessarily comovingly frozen. 
Indeed we show that the domains can expand much faster than the universe, resulting in 
large correlations of the field. 

Here it is important to point out our basic implicit assumption, which concerns the 
damping of the small scale structure of the magnetic field as the correlated domains ex- 
pand. Indeed, although we will not refer to specific damping mechanisms, we assume that 
all the Fourier modes of the magnetic field with wavelength smaller than the dimensions of 
the correlated domains are fully damped and, therefore, the field is coherent inside these do- 
mains. This is equivalent to assuming that the magnetohydrodynamic backreaction, which 
could transfer power to the small scales at later times, is not effective. This assumption 
enables us to avoid the intrinsic non-linearities of the problem (by confining them into the 
damping mechanisms) and to attempt a linear approach without using perturbation theory. 
Of course, a more realistic picture would have to include some transfer of power to the 
smaller than the correlation length scales, especially the ones that are near the scale of the 
correlated domains. In that sense our work can be viewed as the optimum case, i.e. the 
limit of the fastest possible growth of the correlated domains. As such we can still have 
predictive power by setting upper limits to the strength and coherence of the magnetic 
field configuration at any given time. 

In order to describe the evolution of the correlated domains one has to determine the 
peculiar velocity v of equation ([|). This primarily depends on the opacity of the plasma. 

If the plasma is opaque on the scale of a correlated domain, then radiation cannot 
penetrate this scale and is blocked inside the plasma volume. Consequently, the plasma 
is subject to the total magnetic pressure of the magnetic field gradient energy. Therefore, 
this energy dissipates through coherent magnetohydrodynamic oscillations, i.e. Alfven 
waves. This is evident, since the coherent, bulk velocity of the plasma motion is driven by 
the magnetic field so that pv"^ ~ B^. Thus, there is equipartition of energy between the 

^Note that the plasma does not have to be carried from one domain to another or get somehow mixed. 
Also, conservation of flux is not violated with the field's rearrangements, since the field always remains 
frozen into the plasma, which is carried along. 
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coherent motion of the plasma and the magnetic filedQ. Consequently, in this case, the 
peculiar velocity of the magnetic field reorientation, is the well known Alfven velocity [§ ,0 

VA^^ (2) 

where B^d is the magnitude of the magnetic field inside a correlated domain and p is the 
total energy density of the universe, since, before teq, matter and radiation are strongly 
coupled.^ 

Now, if the plasma is not opaque over the scale ^ of a correlated domain, then ra- 
diation can penetrate this scale and carry away momentum, extracted from the plasma 
through Thomson scattering of the photons. This subtraction of momentum is equivalent 
to an effective drag force, F ^ parVTUe [^. Balancing this force with the magnetic force 
determines the "Thomson" velocity over the scale ^, 

VT^J^ (3) 

where va is the Alfven velocity, is the electron number density and is the Thomson 
cross-section. 

Hence, for a non-opaque plasma the peculiar velocity of the plasma reorientation is 
given by @], 

V = min(f A, Vt) (4) 

In order to explore the behaviour of the opaqueness of the plasma, we need to com- 
pare the mean free path of the photon ~ (neCXT)"^ to the scale ^ of the correlated 
domains. For realistic models, the correlated domains remain opaque at least until the 
epoch tank ~ 0.1 sec of electron-positron annihilation (T ~ IMeV). The reason for this 
can be easily understood by calculating It before and after pair annihilation. 

For T > 1 MeV, instead of the usual Thomson cross-section ctt, we have the Klein- 
Nishina cross-section 0, 

3 ,me, r, 2T 1 ^ „ ,GeV . , r ^ i ^ 9 /^n 
a,,^-ar(^)lln- + -]=. 2.7 (— ) ) GeV- (6) 

where rrie ~ O.SGe^is the electron mass and ax — 6.65 x lO^^^cm^ ~ 1707.8 Ge^^^. The 
electron number density is given by jlOi , 



ne^l^-^9eT' (6) 

4 7T^ 



^The coherent plasma motion should not be confused with the thermal motion where Vfh — \JT jm. 
^Unless explicitly specified, natural units are being used {Ti ^ c = 1). 

*This coupling implies that any reorientation of the momentum of matter has to drag radiation along 
with it. This increases the inertia of the plasma, that balances the magnetic pressure. 
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where ^(3) ~ 1.20206 and Qe = 4: are the internal degrees of freedom of electrons and 
positrons. 



From (Kl) and (M) we find, 



It - ^^^7^ for T > 1 MeV (7) 



2^2 

which at annihilation gives, iTiianh) ~ lQ^GeV~^ . 

After annihilation the electron number density is given by ||10||, 



ne ~ 6 X 10^^°n^ ~ 1.44 x lO'^^T^ (8) 



where is the photon number density given by. 



n,:.^g,T^ (9) 



where g-y = 2 are internal degrees of freedom of the photon. 
With the usual value for we obtain. 

It ~ ^3 for T < 1 MeV (10) 

At annihilation the above gives, lT{tanh) ~ 10^^GeV~^ . 

Hence, the mean free path of the photon at the time of pair annihilation is enlarged by 
a factor of 10^°! As a result. It is very likely to become larger than ^ after tank- If this is 
so, the Thomson dragging effect has to be taken into account and the peculiar velocity of 
the plasma reorientation is given by (^). 

In order to calculate the peculiar velocity it is necessary to compute the Alfven velocity, 
which requires the knowledge of the magnetic field value Bed inside a correlated domain. 
To estimate that we assume that the magnetic flux, on scales very much larger than the 
sizes of the correlated domains, is conserved, as implied by the frozen-in condition. 

Consider a closed curve C in space, of lengthscale L > ^, encircling an area A. Con- 
servation of flux suggests that the flux averaged mean magnetic field inside A scales as 
a~^. This implies that for the field inside a correlated domain we have, Bcd{L/C,)~^ oc a~^. 
Since C follows the universe expansion L oc a, with a oc t^^^. Thus, for the radiation era, 
we obtain, 

B,,t'/^^ = K^B,a=jl^ (11) 

where if is a constant to be evaluated at any convenient time. We will show that the 
correlation length grows at least as fast as the universe expands. This implies that the 
magnetic field inside a correlated domain dilutes at least as rapidly as for the radiation 
era. 

Subsituting the above into (0) we find. 
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K 

VAr-lO (12) 

rripi 4 

where we have also used the well known relation, 

t^O.Sg;'/^^) (13) 

where rupi ~ 1.22 x lO^^GeV" is the Planck mass and (7* is the number of particle degrees 
of freedom, which, for T < 1 MeV is 3.36 |g. Q 

Solving the evolution equation (|l|) with a oc t^/^ in the case that v = va gives, 

at? = il)e+'^VA{t)mtii~J^) (m) 



where is the correlation length of the field at the time tj. The first term of (14) is due 
to the Hubble expansion, whereas the second term is due to the peculiar velocity. 

In the case of f = vt, for t > tank, using (^, (|^) and the usual value of ax, (0) gives. 



1 

where 



vt = D— (15) 



D ~ IQ-^'^K^GeV-^/^ (16) 
Using (|TB|), the evolution equation (^ gives. 

The evolution of the correlation length of the magnetic field configuration is described 
initially by the Alfven expansion equation (14) until the moment when ^ ~ /y. From then 
on the growth of ^ continues according either to (|T^) or to (|T^, depending on the relative 
magnitudes of the velocities va and vt- Using the above, the scale ^eq of the correlated 
domains at t^q can be estimated. With a suitable averaging procedure, this this can be 
used to calculate the rms magnetic field over the protogalactic comoving scale at the time 
when structure formation begins. 

3. An important issue, which should be considered, is the diffusion length of the freezing 
of the field. Indeed, the assumption that the field is frozen into the plasma corresponds to 
neglecting the diffusion term of the magnetohydrodynamical induction equation [p!l[] , 

_ = V X (v X B) + a-V2B (18) 

where v is the plasma velocity and a is the conductivity. In the limit of infinite conductivity 
the diffusion term of (|TB|) vanishes and the field is frozen into the plasma on all scales. 



^In natural units G = m^j 



-2 
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However, if a is finite then spatial variations of the magnetic field of lengthscale / will 
decay in a diffusion time, r ~ at^ W^- Thus, the field at a given time t can be considered 
frozen into the plasma only over the diffusion scale, 

~ ^ (19) 

If Id > ^, the magnetic field configuration is expected, in less than a Hubble time, to 
become smooth on scales smaller than 1^. Thus, it is more realistic to consider a field 
configuration with coherence length /^(t) and magnitude of the coherent magnetic field 
Bed, where Bed = Bcd/N is the flux-averaged initial magnetic field over N = 1^/^ number 
of domains. 

An estimate of the plasma conductivity is necessary to determine the diffusion length. 
The current density in the plasma is given by, J = nev, where n is the number density 
of the charged particles. The velocity v acquired by the particles due to the electic field 
E, can be estimated as v ~ eEiTc/m, where m is the particle mass and Tc = Imfp/v is the 
timescale of collisions. Since the mean free path of the particles is given by, Imfp — ^/na^ 
the current density is, J ~ e^E/mfCTc, where dc is the collision cross-section of the plasma 
particles. Comparing with Ohm's law gives for the conductivity |]1T|, [0, 



a ~ —- = (20) 
mvac 47ri/c 

where 

= (^)l/2 is the plasma frequency and z/^ = nvac is the frequency of collisions. 
The collision cross-section is given by the Coulomb formula [0, 

^c^l^lnA (21) 

where In A ~ \n{e~^^T^/n) is the Coulomb logarithm. Thus, the behaviour of the conduc- 
tivity depends crucially on the temperature. 

For low temperatures, T < rrie — 1 MeV (i.e. after tank), the velocity of the electrons 
is, V ~ ^T/nie- Thus, from (|20|) and (^) the conductivity is given by. 




For high temperatures, T ^ m,| (H) suggests that, InA ~ 1. Also, the mass of the 
plasma particles is dominated by thermal corrections, i.e. m ~ T, and w ~ 1. Consequently, 
in this case, (pOD and (0) give for the conductivity, 

a~c^ (23) 

^At not extremely high temperatures the roles of the electrons and the protons may be reversed jl^ 
which would imply that m = rup [nip being the proton mass) instead of me- This would decrease the 
conductivity by a factor of (me/mp)^/^ ~ 0.02. At even higher temperatures T 3> ?Tip, thermal corrections 
become dominant. 
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where c is a numerical factor of order unity.[] 

Using the above resuhs we can estimate the diffusion length. Indeed, from ([l9|) , 
and p3|) we obtain, 

r lO^GeV^/^T-^/^ T > 1 MeV 

k ~ (24) 
[ 10«G'e\/3/42-7/4 J, ^ ^ 

An important point to stress is that the diffusion length is also increasing with time. 
If Zrf > ^ then the size of the correlated domains is actually determined by the diffusion 
length and it is the growth of the later that drives the evolution of the magnetic field 
configuration. 

Another lengthscale necessary to consider is the magnetic Jeans length, given by 



(see also [|T§] and |T3)- This lengthscale is a measure of the dissipation of the field after teq. 
Indeed, in order for the field to have any astrophysical implications in structure formation, 
it needs to be coherent over, 

~ (25) 

(i.e. ^eq > -^s^)- If this is not the case then field oscillations dissipate the energy of the field 
and lead to its decay. If the field is coherent over the above lengthscale then it can have 
an accumulative effect and result in density instabilities, which can, by themselves, lead to 
structure and galaxy formation p!7[] . 

4. Following the above analysis, the strength and coherence of the magnetic field 
configuration at any stage of its evolution^, can be calculated if the intitial values of the 
field and the correlation length are given. 

Being interested in the galactic magnetic fields, we would attempt to calculate the rms 
field on the scale of a protogalaxy ~ 1 Mpc at the time teg, when structure formation begins. 
To do so we need to employ a suitable averaging procedure. 

Choosing the magnetic field as a stochastic variable, Enqvist and Olesen have shown 
that the root mean square value of the field would behave as, 

1 



B,ms = ^{B^) = -7^B,^ (26) 



where is the number of correlation lengthscales over which the field is averaged. At teq 
we have, 

leq 

N=-^ (27) 

^As it is shown in the main contribution to the early universe conductivity is from leptonic inter- 
actions, c is found to be a slowly increasing function of temperature (since at low temperatures most of 
the leptonic species heve underwent a pair-annihilation period). It is shown that c ranges from 0.07 at 
T ~ 100 AleV to 0.6 at T ~ 100 GfeF. For even higher temperatures c approaches the textbook estimate 
c ~ 1.3 |l^ for relativistic electron scattering off heavy ions. 

^though not after t^g 
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where L'g' is the protogalactic scale at teg, for which we find, L'g' ~ {teq/tp^)'^l'^L'PJ ~ lOpc, 
where tp^ ~ lO^^sec is the present time and L^J = 1 Mpc. 

At this point it should be mentioned that, in the above treatment, the rms value of 
the field has been computed as a line average. The averaging procedure could have an 
important effect on the results and has to be considered carefully. One argument in favour 
of line-averaging is that the observed galactic magnetic fields have been measured using 
the Faraday rotation of light spectra, which is also a line (line of sight) computation. If we 
assume that the ratio of the galactic dynamo seed field and the currently observed galactic 
field is independent of the averaging procedure then this would suggest line-averaging is 
required for the computation of the primordial rms field. However, the nonlinearity of 
the dynamo process as well as the rather poor knowledge we have for galaxy formation 
make such an assumption non-trivial. In any case, apart from the above, there seem to 
be no other argument in favour of a particular averaging procedure. Therefore, using 
line-averaging could be the safest choice. 

Here it is also important to point out that line averaging just gives an estimate of the 
rms field and does not correspond to any physical process. Thus, it should not be confused 
with flux averaging which does correspond to a physical process, that of the field untangling, 
and is so in order to preserve flux conservation on scales larger than the diffusion length. 

5. If the magnetic fleld is produced at very early times e.g. during an inflationary era, 
electroweak uniflcation needs to be taken into account. If the original fleld is created before 
the electroweak transition then, assuming that it becomes "frozen" into the electroweak 
plasma is non-trivial. Indeed, during the electroweak era, since the electroweak symmetry 
group SU{2) X U{1)y is unbroken, there are four apparent "magnetic" flelds, three of 
which are non-Abelian. It would be more precise, then, to refer only to the Abelian 
(hypercharge) part of the magnetic fleld, which satisfles the same magnetohydrodynamical 
equations as the Maxwell fleld of electromagnetism. The non-Abelian part of the fleld may 
not influence the motion of the plasma due to the existence of a temperature dependent 
magnetic mass, ms ~ 0.28g^T |]18|, which could screen the fleld. Then, the motion of the 



plasma is determined primarily by the Abelian fleld and can reach a selfconsistent pattern, 
which will "lock" onto the fleld in the same way as in electromagnetism. 

The condition for this screening to be effective is obtained by comparing the screening 
length r5 ~ m^^ of the non-Abelian magnetic flelds with the Larmor radius of the plasma 
motion ~ where m ~ y/aT (a = (7^/47r) is the temperature induced physical mass 
of the plasma particles, g is the gauge coupling (charge) and v is the plasma particle 
velocity. Assuming thermal velocity distribution, i.e. mv"^ ~ T, suggests. 



i?= — -10-2— (2J 



— ~ 10-2 — 
rs Bed 

If i? > 1 then the restriction to the Abelian (Hypercharge) part of the magnetic fleld 
is justifled. This restriction would not cause any signiflcant change to the flnal magnitude 
of the "electromagnetic" magnetic fleld, since, at the electroweak transition, the hyper- 
charge symmetry projects onto the photon through the Weinberg angle, cos6'vk ~ 0.88. If, 
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however, R < 1 then the non-Abehan fields do affect the plasma motion and should be 
taken into account. Since, T oc and B^d falls at least as rapid as a~^, R is in general 
an increasing function of time. Thus, the constraint has to be evaluated at the time of 
creation of the magnetic field configuration. 

6. In determining the rms value of the field at teq one has also to take into account a 
number of constraints regarding its strength. An obvious requirement is that the magnetic 
field should not dominate the energy density of the universe. The expansion of the universe 
dilutes the energy density of the magnetic field, pb = B^JSti, inside a correlated domain 
more effectively than the radiation density, which scales as a~^. Therefore, it is sufficient 
to ensure that psit) is less than the energy density p{t) of radiation at the time of the 
formation of the magnetic field configuration. Thus, the constraint reads, 

PBitd<p{U)^Bl,<^^ (29) 

Another constraint comes from nucleosynthesis. This has been studied in detail by 
Cheng et al. 0. They conclude that, at tnuc ~ 1 sec, the magnetic field should not be 
stronger than. 



Qnuc < IQllGa-uss 



(30) 



on a scale larger than ~ lO^cm. A more recent treatment by Kernan et al. relaxes 
the bound by about an order of magnitude, 5"-'"'^ < e~^(T""'^)^ ~ lO^^Gawss, where is 
the neutrino temperature and e is the electric charge. This bound is valid over all scales. 
Similar results are also reached by Grasso and Rubinstein . 

Finally, an additional consideration would be the lower bound due to the galactic 
dynamo requirements, which, as already explained, demands a field of magnitude. 



> iQ-^^Gauss (31) 

over the comoving scale of 1 Mpc. 

Our field evolution mechanism results in a more coherent magnetic field at teq than 
previously considered. As a result the rms field produced by the mechanisms in could 
also be of greater strength. In some cases, the achieved rms magnetic field could be strong 
enough to dispense with the galactic dynamo. 

At this point we should mention that there is some scepticism regarding the galactic 
dynamo [^, mostly due to the fact that, until now, there is no consistent dynamo model 
pT| . If there is no dynamo mechanism, the growth of the galactic magnetic field is only due 



to winding and, therefore, is linear (in contrast to exponential). Some authors believe that 
this linear amplification of the field would be enough to account for the currently observed 
galactic magnetic field, considering also the additional enchancement by line dragging 
during the gravitational collapse. It has been argued that, if this is the case, a seed field of 
the order of, B^'' ~ lQ~^Gauss would be sufficient. The gravitational collapse would also 
sweep the intergalactic field lines into the galaxies leaving a relatively weak intergalactic 
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field, in agreement with observations, which find that the intergalactic field is less than 
10~^Gauss. 

7. To demonstrate the efficiency of the mechanism, in terms of magnitude and co- 
herency, we will apply it to two toy-models of primordial magnetic field generation. For 
each one of these models the initial correlation length and the time of formation of the 
initial primordial magnetic field configuration is specified in such a way that it could cor- 
respond to a realistic situation. The initial magnitude of the generated magnetic field, 
though, is treated as a free parameter. In both cases we generate the field at a phase 
transition, when the universe is out of thermal equilibrium and, thus, the creation of such 
a field is acceptable p3[|. 



CASE 1. At the electroweak phase transition. 

In general it is thought that the electroweak transition is weakly first order, that is it 
occurs through bubble nucleation. A natural choice for the correlation length of a magnetic 
field, that is created at the transition, would be the bubbles' size at the time of their 
collision. Using typical values for the parameters, we find, ~ (f^/F)^/^ ~ 10^'-'GeV~^ 
10'^ H-^ H, where ~ 0.1 is the velocity of the bubble walls, F ~ IQ-^^GeV^ is the 



bubble nucleation rate per unit time per unit volume EM and H^^ ~ t^w ~ lO^^GeV ^ is 



the horizon size at the time of the transition. We assume that, by some mechanism, an 
initial magnetic field, B^^ ~ lO^Gauss is generated at the transition, where Z is a free 
parameter. Using the above initial conditions we can explore the evolution of such a field. 
From flTTI) we find that, 



K ~ lO^-^GeV'^^ (32) 
Inserting the above into ([12|) we find, 

, +1/2 

-1/2 ^ 



VA ~ 10^-22^^^-1/2 t_ (33) 



Using this in (0) we can estimate the correlation length ^anh at the time tank of pair 
annihilation, 

r lO^^GeV-^ Z <17 

Uh ~ (34) 
[ 10^/2+7^6^-1 z>17 

The above suggests that the Alfven expansion dominates only for Z > 17. Otherwise 
it is the Hubble term that determines the evolution of the correlation length. Using ( p^ 
we can compute the diffusion length at tank- We ffnd that ^anh ^ ^3"^ ~ lO^'^GeV^"^, i.e. 
the correlation length is always larger than the diffusion length at that time. 

From (|12D and (|15]) it can be easily verified that, for all Z, VA{tanh) > vritanh)- Thus, 
from the time of pair annihilation the correlation length evolves according to (|13). Follow- 
ing its evolution likewise it can be shown that the Thomson expansion cannot compete with 
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the Hubble one. Thus, the evolution of the correlation length is driven by the expansion 
of the universe. At teq we find, 

IQ'^^GeV-^ Z <17 

ieq ~ { (35) 

10^/2+i3GeV^-i Z > 17 



The diffusion length at teq is found by p^) to be, fj' ~ lO^^GeV^^. Thus, the above 
correlation length is larger than the diffusion length only ii Z > 20. Therefore, the actual 
size of the correlated domains at teq is given by, 

ieq ~ \ (36) 



With the use of the above and (|Tl]) we find, 

r IQ^-'^^Gauss Z <20 



Z > 20 



(37) 



Also, from (1271) we have 



Z < 20 

N ~ { (38) 

Z > 20 




Using the above (|26|) gives. 



->eq 
rms 



10^^/^-^^ Gauss Z>20 



(39) 



By enforcing the various constraints we can specify the acceptable range of the final 
rms magnetic field strength and the corresponding necessary initial field. 

The requirements of the galactic dynamo suggest that Z > 10. The energy density 
constraint demands that Z <2Q and the nucleosynthesis constraint requires Z < 28. Thus, 
the acceptable range of values for the magnetic field is. 



10 <Z< 26 

(40) 

IQ-^'^Gauss < B^l^ < lO'^Gauss 

Checking with the magnetic Jeans length we also find that the magnetic field could 
influence the structure formation process only ii Z < 13. Finally, from ( ^OD it is apparent 
that for some parameter space the galactic dynamo is not even required since an rms field 
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of strength up to a jiGauss at teq can be achieved. In fact, we can dispense with the 
galactic dynamo for Z > 21. 

CASE 2: At grand unification. 

According to most scenaria for the breaking of grand unification, this occurs at a 
temperature Tqut ~ IQ^^GeV . Since not many details are specified for the nature of 
this phase transition the safest choice of a typical lengthscale is the horizon itself. Thus, 
for the initial correlation length we choose, ~ ~ lQ~^^GeV~^ . Similarly with the 
electroweak treatment, we assume that, at the GUT transition^ a magnetic field of strength, 
~ 10^ Gauss is generated, where W is treated as a free parameter. 

In the same way as in the electroweak case we have, 

K = IQ^-^^GeV^I^ (41) 
with the use of which the correlation length at pair annihilation is, 

r /^'^ ~ l{}^^GeV-^ W <b2 

Uh ~ (42) 
[ 10W2-i3^gy-i W>52 

Again, after pair annihilation it can be shown that VAitanh) > VTitanh)- Thus the 
evolution of the correlation length continues with Thomson expansion. This time, however, 
the Thomson velocity is high enough to dominate the Hubble expansion rate. The resulting 
corrlation lengthscale at teq is, 

r rj' ~ WGeV'^ W <hA 

ieq - (43) 
[ 10^/2-4(^g^-l W>U 

Using the above we can calculate -B^^ and in a similar way as in the electroweak 
case. Then it can immediately be shown that, 

r IQ^-'^^Gauss < 54 

BZs - (44) 
[ l^)^^ I ^-^^ Gauss > 54 

By employing the various constraints on the above result we can determine the accept- 
able range of the magnetic field rms values. The galactic dynamo requirements suggest, 
FT > 49, the energy density constraint demands, W < 54, the nucleosynthesis constraint 
requires, W <QQ and the non-Abelian constraint sets, W < 50. Thus, the acceptable 
range is very small, 49 < < 50 and corresponds to an rms magnetic field of order, 
^rms ~ lO^'^^^^^^^Gawss, which could just about seed the galactic dynamo mechanism. Of 
course, one can imagnine that, if the non-Abelian constraint is violated this should not 
necessarily mean that there will not be any surviving magnetic field. For instance, it is 



'GUT here stands for Grand Unified Theory. 
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highly probable that, if the non-Abelian fields do influence the plasma motion they could 
perturb it in such a way that the magnetic fleld strength would diminish enough for the 
Larmor radious to increase over the non-Abelian screening length (with some of the mag- 
netic energy being thermalized into the plasma). Thus, in this way stronger magnetic 
fields could survive the non-Abelian era. Moreover, such a mechanism could "cut-down" 
the initial strength of a primordial field and reduce it enough never to challenge the nucle- 
osynthesis and energy density constraints. In that sense the non-Abelian stage of evolution 
of the field could expand the upper bound for W instead of reducing it to as low as 50. 
Still, it is difficult to imagnine a way of producing a stronger than 10~^^Gauss field, and 
in that sense the GUT-transition provides a narrow window towards a galactic seed field. 
By checking with the magnetic Jeans length, though, it can be shown that such a field 
would also infiuence the structure formation process and, therefore, it could have additional 
astrophysical effects. 

The above was just a toy-model analysis of our mechanism. The mechanism can be 
applied to a variety of models. A more detailed and complicated example of such an 
application is given in , where a primordial magnetic field is created at the breaking of 
grand unification during infiation, in such a way that the narrow window of the GUT-case 
is achieved in a natural and realistic way. 

8. In conclusion, we have analysed the evolution of primordial magnetic fields and 
shown that, when the effects of the surrounding plasma are taken into account, the corre- 
lation length of the field configuration grows faster than the Hubble expansion. This results 
into a more coherent magnetic field than previously thought. There is some similarity with 
the recent work of Brandenburg et. al. who describe the plasma with relativistic 



magnetohydrodynamic equations and turbulent behaviour. However, they are unable to 
solve their equations analytically, and, instead, they attempt a numerical analysis in 2+1 
dimensions. This also results in a faster growth of the coherence length. Still, they do not 
take into account all the effects of the plasma on the magnetic field configuration, such as, 
for example, the Thomson scattering effect. 

We would like to thank K. Bajers, B. Carter, N. Turok, T. Vachaspati and especially 
M.J. Rees for dicussions. This work was partly supported by PPARC, the Greek State 
Scholarships Foundation (I.K.Y.) and the E.U. under the HCM program (CHRX-CT94- 
0423). 
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